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We revisit the definition of the probability current for the Schrodinger equation. First, we prove 
that the Dirac probabihty currents of stationary wave functions of the hydrogen atom and of the 
isotrop harmonic oscillator are not nil and correspond to a circular rotation of the probability. 
Then, we recall how it is necessary to add to classical Pauli and Schrodinger currents, an additional 
spin-dependant current, the Gordan current. 

Consequently, we get a circular probability current in the Schrodinger approximation for the 
hydrogen atom and the isotrop harmonic oscillator. 



I. INTRODUCTION 



In this letter we revisit the definition of the probability current for the Schrodinger equation. We discuss the 
probability current of the hydrogen atom and the isotrop harmonic oscillator for which the classical results can be 
expressed in the form: "the probability current of the wave eigenfunctions ipnim of Schrodinger with m = is nil", 
cf. for example Indeed, as the function V'nio = Rniir)Yio{9,ip) is real, the classical definition of the Schrodinger 
probability current leads to 

Ji = :^(W*-^*VV;) =0. (1) 

First, we prove that for the Dirac equation, the probability current of stationary wave functions of the hydrogen 
atom and of the isotrop harmonic oscillator are not nil and correspond to a circular rotation of the probability. 

Then, we recall how a good approximation of the Dirac probability current is obtained in the nonrelativist case. It 
is necessary to add in equation ^ the Gordon current 2] which reads 

J2 = -^rot{'^* ail)) 
Im 

in the Pauli approximation |^, Q and 

J2 = — Vp X s (2) 

TO 

in the Schrodinger approximation P|, where s corresponds at a constant spin vector. Consequently, we get 
a circular probability current in the Schrodinger approximation for the hydrogen atom and the isotrop harmonic 
oscillator. 



II. THE DIRAC PROBABILITY CURRENT OF THE HYDROGEN ATOM 

In the presence of electromagnetic couplings, the Dirac equation reads 

-^7"^ + |7'^^.^ + ^V' = (3) 
where the 7'^ (/i = 0, 1, 2, 3) are the Dirac matrices. 
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The probability density p = tp'^'tp and the three components of the probabihty current J: 
satisiy the continuity equation 



dp 
dt 



V.J 0. 



(4) 



(5) 



In the case of an electron in the hydrogen atom with a spherically symmetric potential V = and a potential 

vector A = 0, the Dirac eigenvalues equation is then 

E cV Tnc 
M-j- + 7-)V'(i-) - iT.ViPir) + -j-^iv) = 0. 

The eigenfunctions ^n.i,j.m{'^) depend on the quantum numbers n,l,j and m with n integer (n > 1), I integer 
(0 < I < n — I), j ~ ^ if I — and j — I + m if I > 1, m — ±i. Then il'n,i,j,m{'^) is given y| by 

with /' = J + i if / = j — i and I' = j — \ ii I = j + \ , and where /(r) and g{r) are real functions of r and where 



21 + 1 ''l,m-^ 



21 + 1 ^l,m+^ 



We obtain four classes of wave eigenfunctions: m = ^ or — i, / = or / > 1. For ni = ^ 



( 



'^n,l,],m=+iir,e,ip) 



fir) 
ig{r) 



\ 



a{e) 

b{9) exp(j(^) 
d{9) exp{i(p) J J 



(6) 



with for / = 0: a{e) = Yqo, b{9) = 0, c{e) ^ -Yoocose, d{9) = -yoosin6' and, for / > 0: a{9) = y^^F/^o, 



^(^) = \/ STTT^M cxp(-«</3), c{9) = y/ 5j^i1+i,o and = cxp(-z<y9), and for m = -i 



f d{9) cxp(— 
I ~c{9) 



ig{r) 



(7) 



with for I ^ 0: 0(6*) = 0, b{9) = Foo, c{9) = -roosin6', d(6l) = roocos6l and, for / > 0: a{9) = -^^^Y;,-! exp(i(^), 
K^) = c[9) = -yi^l^_i__iexp(i(^) and = ,J^Yi^,^o- 



For a wave function 



and the Dirac matrices 



/ 

-/ 



and 



a' 




the probability density is equal to: p = XiXi +X2X2, and the probability current is given by: J' = c{Xi'^^X2+X2^^Xi)- 
For the functions ® and Q, we find : p = p{r)[a^{9) + 9{9)]+g'^{r)[c^{9)+d'^{9)], = -2cf{r)g{r)[ad~bc]smip, 
,P = —2cf{r)g{r)[ad — be] cos 93, 0, and as probability current velocity: 



J 

V = — = 



2cf{r)g{r)[ad-bc] 



p p{r)[a^ +b^]+ g^{r)[c^ +d^] 
which corresponds to circular probability currents. 



(8) 
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We have the same resuh ^ for all the Dirac equations with a spherically symetric potential; in particular for the 
isotrop harmonic oscillator. 

The Dirac eigenf unctions are not stationary. All the electron wave eigenfunctions of the hydrogen atom correspond 
to circular probability currents. It is not the case for the Schrodinger equation with the current given by (Q, since 
we have then stationary solutions for the eigenfunctions 'tpnio- In the fondamental state Isi with the energy E'^ i = 

meC^VT^, we have: /(r) = ao{^)^^-^e'^ , g{r) = -a^ ^--^ {l-)^^~^e'^> , where = ^ = 0.53 A 

is the Bohr radius. From ®, we deduce p{r,e,Lp) = p{r) ^ 2a2i^:^^^( J-)2(VT^-i)e"^^ , and 

V — ac sin 0Uy. (9) 

From O and ©, we deduce v = —acsvaOu^. There is a big difference between and the current nil given 
by 1^ with the Schrodinger equation. The lesson of this example is that, although Schrodinger equation is a good 
approximation of the Dirac equation, the Schrodinger probability current Q is not the good approximation of the 
Dirac probability current This conclusion appears clearly, as we recall now 0]) P]i in the approximations to 
transform first the bispinor of the Dirac equation to the spinor of the Pauli equation, then the Pauli spinor to the 
wave function of the Schrodinger equation. 

III. THE PAULI AND SCHRODINGER PROBABILITY CURRENTS 

After the change of variable ip = ip's'T^"^'^ *, the Dirac equation J^J for the bispinor tp — (xijX2) can be written 

[ih— ~ eV]xi = ca{-ih\/ - -A)x2 (10) 



d e 
[ih— - eV + 2mc^]x2 = ca[~ihV - -A)xi. (11) 

Small electron velocities allow us to neglect the two first terms of the left part of equation Hll|) and to eliminate X2 
in H10() . So we obtain the Pauli equation for the spinor xi 

where H = V x A is the magnetic field. 

As X2 ^ Xii the Dirac density p = ijj'*ip = XiXi + X2X2 is equal in first approximation to the density p = XiXi of 
the Pauli equation. Then the Dirac current density can be written 

J = il^iXi^xl - xl^Xi) ~ —Axixl + t/^V X ixl<7Xi) (12) 
Zm mc Zm 

Thus, to obtain in Pauli a good approximation of the Dirac probability current, it is necessary to add the Gordan 
current 

to the classical Pauli current 

h e 

Ji = iT^(xiVxJ - XiVxi) AxiXi- 

zm mc 

Then, the current J = Ji + J2 verifies, as the current Ji, the continuity equation(|SJ. 

To obtain the Schrodinger equation from Pauli equation, consider the case where we have no magnetic field Q and 
where the system is in a spin eigenstate 

Xi(r,i) = ip{r,t)x 

with a constant spinor x such as x*X = 1- Then the function ip{r, t) verifies the Schrodinger equation 

.^.-l-A^ + ey^. (13) 
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However, in this case, the current H12(l does not reduce [Hj to classical Schrodinger probability current. Indeed, 
writing ip{r, t) — yfpe^Ti , we have 



with 



J = — VS'+ — V/9 X s 
m m 



Indeed, to obtain a good approximation of the Dirac probability current for the Schrodinger equation, it is necessary 
to add the spin-dependent current 



to the classical Schrodinger current 



— Vp X s 

771 



Ji = — VS'. 

777 



Therefore the probability current depends on a constant spin vector s. The add of a constant spin vector to the 
Schrodinger equation is not new, and has been used before by Landau J^. 

This current J2 gives also a simple explanation of the value 2 of the gyromagnetic ratio of the Dirac electron; 
Indeed, this contribution to the particle orbital angular momentum is 

L2 = r X mJ2 = r X Vp x s 

and this mean angular momentum is then 

< L2 >= j l.-id^r = J{r ■ s)Vptfr - J (r ■ \7p)sd^r = -s + 3s = 2s. 

We can verify the precedent theoritical approach with the computation of the current H13|) of the hydrogen atom 
wave eigenfunctions in the approximation of Schrodinger. 

For the wave functions il'rdm{'>',0,'^) = i?„;(7')P;„j(0) exp(i777(/3) exp(— i-^^), the classical Schrodinger current is 

p mh 

Ji = 

777p r sm U 



with p = Rni{'r)Pf„^{0). The current is nil for the wave functions with 777 = 0, in particular for the states ns. If 777 > 0, 

fkands = -fl 

s — sgn{m) (cos6'Ur — sin^u^^) 



we take as spin vector s = §k and s = — $k if 777 < 0. Then, this vector can be written 



and the additional spin-dependent current (Gordon current) J2 is then equal to 

h 

J2 = -sgn{m) — (i?„ii?;,;P;2^j sin 6* + ^Pi™P/,„ cos 6l)u 

777e r 

Then, the probability current velocity is 



m 



V = - = ^[~sg77 ( 777)(^-^ Sine + cos9) + -^]u^ 

p me tiniyr) rFim[0) rsmO 

which corresponds to circular probability currents exactly as for the Dirac equation (jH)) . In the case of the fondamental 
state Is, we find again exactly the velocity of the Dirac probability current 

Vis = acsinSUip. 

For the states 2s, 2po, 2pi and 2p_i, we obtain for the current velocity 

•^c,, 1 . . „ 

V2s = YT"^)sm0u^, 

2ro 

C^C QIC 

V2po = V2pi = — sin^u^ 0770? V2p„i = -— sinSu^. 
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IV. CONCLUSION 



As we have shown, although the Schrodinger equation is a good approximation of the Dirac equation, the Schrodinger 
probabihty current Q is not the good approximation of the Dirac probabihty current The main conchision is 
that it is necessary to add to the classical Schrodinger current the spin-dependent current 

J2 = — Vp X s 

m 

which corresponds to a constant spin vector s. 

In the case of the isotrop harmonic oscillator with ^(r) = imtj^r^, the fondamental wave function Is 

, , .muj .3 muj 9. , SiLut. 

V'0,0.0 r = i^)' ^''P "^"^ '^'^P ^ 

TTh zfi 2 

gives a Schrodinger classical current Ji nil. If we consider an oscillator with a constant spin s = ^h, then the 
equation gives a circular probability current J with the velocity 

Vis = — = —Lor sinOu^. (14) 

It is the same velocity as in the classical case. For the other eigenf unctions of the harmonic oscillator, we find also a 
non nil probability current, but non circular as in H14|) . 
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